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Abstract. 

The existence of the initial value constraints means that specifying initial data 
for the Einstein equations is non-trivial. The standard method of constructing initial 
data in the asymptotically fiat case is to choose an asymptotically flat 3-metric and 
a transverse-tracefree (TT) tensor on it. One can find a conformal transformation 
that maps these data into solutions of the constraints. In particular, the TT tensor 
becomes the extrinsic curvature of the 3-slice. We wish to understand how the physical 
solution changes as the free data is changed. In this paper we investigate an especially 
simple change: we multiply the TT tensor by a large constant. One might assume that 
this corresponds to pumping up the extrinsic curvature in the physical initial data. 
Unexpectedly, we show that, while the conformal factor monotonically increases, the 
physical extrinsic curvature decreases. The increase in the conformal factor however 
means that the physical volume increases in such a way that the ADM mass become 
unboundedly large. In turn, the blow-up of the mass combined with the control we 
have on the extrinsic curvature allows us to show that trapped surfaces, i.e., surfaces 
that are simultaneously future and past trapped, appear in the physical initial data. 
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1. Introduction 

When we consider General Relativity as a dynamical system, or try to construct 
solutions to the Einstein equations numerically, we draw on our understanding of and 
experience with other theories of physics. As in the case of particle mechanics or 
electromagnetism, we can pose initial data, and use the field equations to propagate 
the system into the future (or even into the past if we wish). We need to give a 
Riemannian 3-metric gij, and a symmetric tensor K % i that is the extrinsic curvature of 
the 3-geometry. All regular solutions to the equations can be obtained in this way. 

Just as in electromagnetism, the initial data must satisfy constraints. This means 
that we must first specify 'free' data, and then manipulate these data to map them into 
a solution of the constraints. What happens to the physical solution as one changes 
the free data? In this article we change the free data in a particularly simple way and 
investigate the consequences. 

Since electromagnetism is a massless spin-1 field theory, and gravity is a massless 
spin-2 field theory, it is not surprising that we can see parallels between them. Consider 
the Maxwell free data as a pair of 3- vectors, (A, F). A is the magnetic vector potential, 
and F is the vector field from which we extract the electric field. Thus, the magnetic 
and electric fields can be generated by B = V x A and E = F — V V, where V is a scalar 
function chosen to satisfy V 2 V = ViF\ Electromagnetism as a theory can be expressed 
in terms of (A, E) rather than the standard (B, E). We can see that A is analogous to 
the metric, since both are gauge dependent quantities, while E is the analogue of the 
extrinsic curvature, because E is essentially the time derivative of A, and the extrinsic 
curvature is essentially the time derivative of the metric. Thus F is the analogue of the 
free data from which we extract the extrinsic curvature. It is clear that if we multiply 
F by any constant, then E will then be multiplied by the same amount. As a result, 
the electromagnetic energy density, B 2 + E 2 , increases without limit, and then the total 
energy will also blow up. In this article we want to multiply the free data, which specifies 
the extrinsic curvature, by a large constant in order to see what happens. We will again 
show that the total gravitational energy increases without limit. 

Initial data for the Einstein equations consists of two parts: (1) a 3-dimensional 
spacelike slice with a Riemannian metric, g^, and (2) a symmetric 3-tensor K l K The 
spacelike slice is to be regarded as a slice through a 4-dimensional pseudo-Riemannian 
manifold that satisfies the Einstein equations, and K l i is the extrinsic curvature of the 
3-surface embedded in the 4-manifold. This means that K 1 ^ is the Lie derivative of 
the 3-metric along the timelike normal to the slice. These initial data cannot be freely 
specified; rather, they satisfy 4 constraints, which, in vacuum, read 



where is the scalar curvature of the 3-metric, K = g a bK ah is the trace of the 
extrinsic curvature, and V* is the covariant derivative with respect to g^. Eq.([T]) is 



(3) R - K ij Kij + K 2 =0 
V i (K ij -g ij K) = 0, 



(1) 
(2) 
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called the Hamiltonian constraint, and Eq.(J2]) is called the momentum constraint. A 
comprehensive account of the constraints can be found in pQ, especially in Chapter 
VII; the classic article is [2]. If K = 0, the momentum constraint shows that 
is both divergence-free and trace-free. A symmetric 2-tensor is called a TT tensor if 
it is simultaneously divergence-free and trace-free. It turns out that TT tensors are 
conformally covariant in the sense that if F? T is TT with respect to a metric gr^, then 
(fi~ 2 FT T is TT with respect to g^ = 4> A 9ij- 

The standard way of finding sets (gij,K^) that solve the constraints is the so- 
called conformal method which was initiated by Andre Lichnerowicz |3J. The idea is to 
specify a base-metric gij and a TT tensor. The base-metric is conformally transformed 
via cjij = <p 4 gij, choosing the conformal factor <p so as to satisfy simultaneously the 
Hamiltonian constraint and the momentum constraint. 

This article deals specifically with constructing initial data that are asymptotically 
flat. This means that we choose the base metric to be asymptotically flat, and seek 
a conformal factor, 0, such that <fi — y 1 at infinity to maintain asymptotic flatness. 
As discussed above, we will restrict our attention to 'maximal' slices, initial data for 
which K = 0. This is a non-trivial assumption: there exist spacetimes which contain 
no maximal slice, see, for example, |3]. On the other hand, there exist many spacetimes 
that do possess maximal slices. The maximal assumption is useful because it makes it 
easy to specify the 'free' data. 

In Section II, we will see what happens when we multiply the 'free' TT tensor by a 
large constant. This maintains TT-ness of the tensor, and we can always find a suitable 
conformal factor to solve the constraints. We will show that, despite the fact that we 
continue to demand <p — > 1 at infinity, the conformal factor monotonically blows up as 
we increase the TT tensor. Further, we can show that the conformal factor behaves like 
(j) m 1 + E/2r + . . . at infinity, where E is the ADM energy of the solution. In Section III 
we will show that E increases monotonically so that the energy becomes unboundedly 
large. In Section IV we show that trapped surfaces must appear in the physical data 
when the scaling is large. 

This paper is a companion to [5], in which we did the same analysis on a compact 
manifold. In [5] we found that the behaviour of the conformal factor was complicated; 
we got blow-up but not uniform blow-up. That paper therefore used a combination of 
numerical and analytic techniques to describe how the conformal factor changed. We 
assumed we would have similar difficulties in this calculation and did some numerical 
modeling. However, on viewing the numerics it was clear that the the conformal factor 
behaved in a much more regular fashion: we clearly had uniform blow-up. This inspired 
us to go back and look again at the analytical work. We were able to produce exact 
mathematical proofs of everything we wanted, and so we were able to eliminate the 
numerics. This, therefore, is a paper that has passed through a numerical phase and 
emerged finally clS du purely analytic document. 
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2. Scaling the Extrinsic Curvature 

The Hamiltonian constraint is 

(3) i? - K ij Kij + K 2 = 0. (3) 

We assume that we are given a suitably asymptotically flat metric, g^, and a suitably 
asymptotically flat TT tensor, KjLp, with respect to the given metric. We map this 
set of initial data to an asymptotically flat solution of the Hamiltonian constraint by 
using glj = (p^gij, and Kjl T = (p~ w K^ T , and finding an appropriate 0. Since the 
extrinsic curvature stays TT under this transformation, the momentum constraint is 
automatically satisfied. This reduces to solving the Lichnerowicz equation 

V 2 - f + \a^- 7 = (4) 

where A 2 = K^ T K^ T . We impose the boundary condition — y 1 at infinity. 

If we have a maximal solution, i.e., K — 0, then the Hamiltonian constraint 
simplifies slightly to ™R — K^Kij = 0, and an immediate consequence is > 0. 
If we have a solution to the Lichnerowicz equation Eq. (J4]) we will conformally transform 
the base metric into one with non-negative scalar curvature. It turns out that there is an 
obstruction to this. Riemannian 3-manifolds split into three classes, called the Yamabe 
classes [TJ. Only metrics in the positive Yamabe class can be conformally transformed 
into metrics with non-negative scalar curvature. These are the metrics that can be 
conformally mapped into a closed, without boundary, compact manifold with constant 
positive scalar curvature. All we need is that the AF metric be conformally flat at 
least to order 1/r 2 . Obviously, if we want to construct a maximal solution to the 
Hamiltonian constraint via a conformal transformation, the base metric must be in the 
positive Yamabe class. One can show that Eq. (T4|) will have a unique positive solution if 
and only if the metric belongs to the positive Yamabe class [8]. 

We want to change the initial data and see what happens. In particular, we rescale 
the background TT tensor. We change Kj! T to a A Kji T where a is a constant. The new 
K l i is still TT and we continue to be able to solve the Lichnerowicz equation. In this 
article we are interested in the behaviour of the conformal factor as a — > oo. This means 
that we write the Lichnerowicz equation as 

V 2 - f + \a*A 2 ^ = 0. (5) 

We are, of course, particularly interested in the properties of the physical initial data, 
(gij,K^ T ) = ((f) 4 gij, 4>~ 10 a 4 K^ T ), as a becomes large. Since we are in the positive 
Yamabe class, we can always make a preliminary conformal transformation to map the 
metric to one which satisfies R = 0. This conformal transformation depends only on the 
base metric, i.e., it is independent of a. This reduces the equation we wish to analyse 
to 

V 2 <P + la 8 A 2 <f)- 7 = 0. (6) 
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with, of course, — > 1. We know that this equation will have a regular positive solution 
for any finite a. What happens to as a becomes large? This equation would make 
no sense if remained regular and bounded as a — > oo, because we would have a finite 
term equalling a term that becomes unboundedly large. Therefore we expect that 
blows up as a becomes large. More precisely, we expect that blows up linearly with 
a. This is why we choose a A as the rescaling factor. We now introduce a normalized 0, 
= <f)/a. This allows us to rewrite Eq.flHJ) as 

V 2 + -A 2 (f)- 7 = 0, -> 1/a at oo. (7) 
8 

Obviously, depends on a. Let us consider a slightly different equation, 

V 2 ip + -A 2 ^' 7 = 0, HO at oo. (8) 
8 

This equation has a regular positive solution which is unique. This ip, when used 
as a conformal factor, results in a compact without boundary manifold satisfying 
R = ip~ 12 A 2 . We can do this in two stages. First: the Yamabe theorem [7] tells us 
that we can conformally map the asymptotically flat metric to a compact manifold of 
constant positive scalar curvature, R . Second: following this transformation, Eq.© 
becomes 

V 2 - f + \A 2 <p- 7 = 0. (9) 

This has a regular solution on the compact manifold. The solution of Eq.® is just 
the product of the compactifying conformal factor by the solution of Eq.(|H]). Note: this 
requires that A 2 be well behaved on the compact manifold. This is obviously satisfied 
if A 2 falls off like r -12 or faster on the asymptotically flat manifold. 
The key reason for introducing ip is that we will show 

<j> —7- ip as a — > oo. (10) 

First: we can show that <p monotonically decreases as a increases. To see this, 
differentiate Eq.(J7|) by a to give 

V 2 f - 7 -A 2 ^f = 0, f -1/a 2 at oo. (11) 
da 8 da da 

It is clear that d(j)/da is negative at infinity. Let us assume that it is positive somewhere 
in the interior. If it is, then d<fi/da will have a positive maximum. At such a point, both 
of the terms in Eq. dlip are negative, and this cannot be the case. Therefore we have 
d(p/da < 0. This means that monotonically decreases as a increases. 

Second: we can show that — 1/a monotonically increases as a increases. The 
equation that 0' = — 1 ja satisfies is 

V 2 0' + Ia 2 (4>' + 1/a)- 7 = 0, 0'^O at oo. (12) 
8 

Again, we differentiate this with respect to a and get 

V 2 ^-^ 2 (0'+l/a)- 8 ^+-^A 2 (0'+l/a)- 8 = 0, d 4- at oo.(13) 
da 8 da 8cr da 
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Let us assume that d<fr'/da is negative somewhere in the interior. Then there must exist 
a point where it is a negative minimum. At such a point all the three terms in Eq. (1131) 
are positive which is an obvious contradiction. Therefore we have d(f)'/da > 0. This 
means that 0' = — 1/a monotonically increases as a increases. 

Furthermore, we can also show that ip satisfies > ip > 0'. To see this, let us first 
subtract Eq.([8]) from Eq.© to get 

V 2 (4>-ij) + Ia 2 (4>- 7 -ij~ 7 ) =0, 4>-if>-> 1/a at oo. (14) 
8 

Let us assume that (0 — ip) goes negative in the interior. This means that there will be 
a negative minimum. At this point we have V 2 (0 — ip) > and, since < ip we have 
0~ 7 > tp~ 7 . Therefore both terms in Eq. (fT4"]) are positive which cannot be the case. So 
we must have — ip > for all a. 

To show that ip > 0', we need to subtract Eq.( TT2|) from Eq.® to get 

V 2 {^ -lp') + \A 2 [ij- 7 + l/a)- 7 } = 0, (V>-0')->° at oo.(15) 
8 

Let us assume that we have a region with ip < <j>'. In such a region ip~ 7 > 4>'~ 7 > 
(<p' + l/a)~ 7 . In this region ip — <p' will have a local minimum. At that point both terms 
in Eq.( fl5l) will be positive, again a contradiction. This implies ip > (/>'. 

Of course, we have — 0' = 1/a. This gap goes uniformly to zero as a — )■ oo and 
we have > ip > 0'. Therefore and 0' uniformly approach ip, one from above and one 
from below, as a — >■ oo. 



3. Controlling the ADM mass 

We know that the leading term of ip is of the form C/2r and we can easily show 



1 



A 2 tp- 7 dv. (16) 



167T 

This means that C is positive and finite. 

Now will behave asymptotically like D/2r + 1/a with D depending on a. We 
will show that D — > C as a — > oo. We have that 

1 / — ■-■ ' / -72 1 j 1 / 42.7-7, 



d = -s i v -^ s ' = -s y v = lis y ^ _7 *- (i7 » 

We know that is bigger than ip but less than ip + 1/a. Thus 

V> < < </> + 1/a. (18) 

Consequently 

(^ + l/a)- 7 < 0- 7 < V" 7 , (19) 

which implies 

— / A 2 (i) + 1/a)- 7 ^ < — [ A 2 4>- 7 dv < — [ A 2 ip~ 7 dv. (20) 

167T J 167T y 16tt 7 
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This can be immediately rewritten as 

— / A 2 (i) + l/a)~ 7 dv <D <C. (21) 

167T J 

Using the monotone convergence theorem we see that j 1 — J A 2 (ip + l/a)~ 7 dv — > C as 
a — > oo. Therefore we get D — > C. From this we conclude that the coefficient of the 
l/2r part of the 'physical' conformal factor, = a<p, behaves like aD rj aC. In other 
words, it blows up linearly with a. 

The energy of the physical metric is the sum of the monopole part of the 
'background' which is finite and bounded plus the monopole part of the conformal factor 
which is positive and becomes large. Thus the total energy is positive and increases 
monotonically with a. This, by the way, is strong-field positive energy proof. 

The conformal factor which maps to a physical asymptotically flat data set is 
(f), the solution of Eq.(jHJ). We use this to generate a solution of the constraints 
(gij,K% T ) = (0 4 ^-,0- lo a 4 ^ T ). From this we get K^ T K^ T = f 12 a 8 i^ T iff '. Since 
scales with a, this goes pointwise to zero like a -4 . 

The physics here is a bit tricky. We can think of Kji T K^ T as the analogue of 
E 2 in electromagnetism, and so can be thought of as the kinetic energy part of the 
gravitational wave energy density. This shrinks like a~ 4 . However, the physical volume 
increases as a 6 . This is why we see the total energy increasing despite the decrease of 
the energy density. 

Nevertheless, it is clear that the physical extrinsic curvature shrinks as a increases. 
Therefore the physical solution looks more and more like a moment of time symmetry 
data set. Each physical solution will have an ADM energy- momentum. Since we assume 
the extrinsic curvature falls off faster than 1/r 2 , the ADM momentum vanishes, and 
the ADM energy becomes the ADM mass. This will be contained in the 1/r part of 
the physical metric. This will be made up of two parts. One part will be the 'mass' 
contribution of the background metric, gij, which is a fixed number. The other part is the 
coefficient of the l/2r term in the conformal factor, <p. This is aD w aC. Therefore, for 
large a, the ADM mass diverges linearly with a. This is a 'strong-field' positive energy 
proof since we show that the energy is positive for large a without any restriction on 
the background energy. 

Further, we can multiply Eq. (I19j) by a to see that <fi satisfies 

ai/j < <j) < aip + 1. (22) 
Therefore the conformal factor blows up uniformly and linearly with a. 

4. Appearance of trapped surfaces 

Trapped surfaces are defined by negative outgoing null expansions. Consider a closed 
2-surface in an initial data set, with outgoing unit spatial normal, n l . The outgoing null 
expansions are defined by 

H± = k± (g^ - n,jnj)K ij = V l n l ± (g^ - n inj )K lj . (23) 
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Given a moment-of-time-symmetry slice in isotropic coordinates through the 
Schwarzschild spacetime, we know that the sphere defined by r = m/2 is the apparent 
horizon and all surfaces with r < m/2 are trapped. In particular, we can show that the 
mean curvature, k, of the surface defined by r = m/4 equals —2/r = — 8/m. We have 
a sequence of physical metrics, labeled by a. On these metrics we can introduce quasi- 
isotropic coordinates near infinity, and the surfaces labeled by r = m/4 will be included 
in these domains for large enough m. We work out the null expansion of these surfaces. 
This will have a negative leading term, —8/m. Now we show that the corrections all 
fall off faster than 1/r, i.e., faster than 1/m. Thus, for large m, the leading negative 
term dominates and the surfaces are trapped. 

The key term in our analysis is k = Vm\ the mean curvature of the 2-surface in the 
3-space. We can ignore the extrinsic curvature term in this calculation because it falls 
off rapidly both with radius and with a. Therefore to find a trapped surface we need 
only find a 2-surface with negative mean curvature, k. This already shows that we are 
getting trapped surfaces in the language of Penrose. The surfaces will be simultaneously 
future and past trapped. 

To do this we need to understand the behaviour of <fi at large radii. From Eq. (1221) 
we know that we can write <fi as 

^=l + ^ + / 5 (24) 

where / is dominated by the dipole moment, so it falls off like 1/r 2 , and grows 
proportional to a. This means / = 0(a/r 2 ). We also need control of the gradient 
of <fi, i.e., the gradient of /. To do this, let us return to Eq.(j7j) and differentiate it with 
respect to, say, x to get 

V A <p — = S, > at oo. (25) 

ax 8 ax ax 

where S is a source term which comes from differentiating A 2 and the metric. This is a 

nice linear elliptic equation which guarantees that d(f>/dx is well behaved. Since m acfi, 

we can deduce that df jdx = 0(a/r 3 ). 

Let us assume that the base metric is conformally flat outside some region of 

compact support and let us compute the mean curvature of a spherical surface r = r , 

assuming that the physical metric is conformally flat, = 4> A 5ij. The flat space unit 

normal, qi = (x/r,y/r, z/r) is proportional to the normal in the physical space. We 

define \q\ 2 = g^q%qj- The unit normal to the surface is then = qi/\q\- We can now 

compute 

k = = J-diirfrf) = -Ld^g^nJ 
V9 V9 

= (1 + m J 2r + /)6 ^[(l + *D/2r + /) 4 (x/r, y/r, z/r)] 

2(l-aD/2r + f + x%f/2) 

r(l + aD/2r + /) 3 ' 1 ' 
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The mean curvature will be negative if 



1 - aD/2r + f + x'dJ/2 < 0. 



(27) 



We know that 



/ + x%f/2 < C 3 a/r 2 



(28) 



where C3 is a constant independent of a. If we choose r = m/4 = aD/A, we get 



The third term will be less than 1 for a large and so the mean curvature goes negative. 

It is a straightforward exercise to show that if the base metric is not conformally 
flat, but is well behaved near infinity, the correction to k falls of quickly. In the same 
way, as we have indicated above, the correction due to the extrinsic curvature is also 
negligible. 

This article is closely related to those cited in [9]. The space of free data for 
the maximal constraints consists of smooth, Riemannian 3-metrics that have positive 
Yamabe constant together with finite TT tensors. We want to find out what happens as 
we approach the boundary of this space. In the articles cited in [9] we kept everything 
regular as we looked at a sequence of metrics along which the Yamabe constant went 
to zero. In this article we again keep everything regular but let the TT tensor diverge. 
This is a very different part of the boundary of the space of free data. Nevertheless, we 
find very similar behaviour. The conformal factor blows up, the ADM mass diverges to 
+00 and horizons appear. 

Unfortunately, this is not enough to construct a valid positive energy proof because 
one needs to investigate the other parts of the boundary, for example we need to consider 
sequences of metrics which either become 'rough' or cease to be Riemannian. 

5. Appendix: Showing D — > C 

It is 'obvious' that as a — > 00, J A 2 (ip + l/a)~ 7 dv — > C, but we have difficulty 
proving it. Here we will prove something weaker, but this result is really all we need. 
Pick a radius r = ro such that 



where B(ro) is the ball inside r = r$. We know that both ip and tq are independent 
of a. Since tp is positive so it will have a minimum value inside B(ro). We pick cto 
so that l/«o equals this minimum value. Therefore, for all a > «o, we have that, for 
r < r , t/> > I /a. Therefore for r < r we have ip + 1/a < 2ip. In turn we get 



1 - aD/2r + f + x l dif /2 < 1 - 2 + 



(29) 



aD 




(30) 




(31) 
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This means that we can write Eg. (121 p as 

£ < D < C. (32) 

We can adjust this proof so as to replace the 2 s by a number which is as close to unity 
as we wish. Pick a /3 << 1. We move r out until we get 

-L / A^~ 7 dv > (33) 
167T J B(ro) 1 + /3 

We also increase a such that for all a > a we have that, for r < r , /?■?/> > 1/a. 
Therefore for r < r we have if) + 1/a < (1 + /3)i[). Then we get 

(TTW <D<G (34) 

This equation now holds for all a > o.q. In other words, we can make D as close to C 
as we wish by picking a large a. 
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